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ABSTRACT 

The flat spectrum radio quasar 3C 454.3 underwent an extraordinary outburst in 
December 2009 when it became the brightest 7-ray source in the sky for over one week. 
Its daily flux measured with the Fermi Large Area Telescope at photon energies E > 100 
MeV reached -Fioo = 22 it 1 x 10^^ ph cm~^ s~^, representing the highest daily flux of 
any blazar ever recorded in high-energy 7-rays. It again became the brightest source in 
the sky in 2010 April, triggering a pointed-mode observation by Fermi. The correlated 
7-ray temporal and spectral properties during these exceptional events are presented 
and discussed. The main results show flux variability over time scales less than 3 h 
and very mild spectral variability with an indication of gradual hardening preceding 
major flares. No consistent loop pattern emerged in the 7-ray spectral index vs flux 
plane. A minimum Doppler factor of ~ 15 is derived, and the maximum energy of a 
photon from 3C 454.3 is ~ 20 GeV. The spectral break at a few GeV is inconsistent 
with Klein-Nishina softening from power-law electrons scattering Lya line radiation, 
and a break in the underlying electron spectrum in blazar leptonic models is implied. 



Subject headings: Galaxies: active 
observations 



quasars: individual: 3C 454.3 — gamma rays: 



Introduction 



The radio source 3C 454.3 is a well-known flat spectrum radio quasar (FSRQ) at redshift 
z = 0.859. It entered a bright phase starting in 2000, and has shown remarkable activity in 
the past decade. It underwent major outbursts in 2005, reaching an R-b and magnitude of 12.0 



and the largest apparent optical lumi nosity ever recorded from a blazar (IFuhrmann et al. 



Villata et al. 



2006 



2007; Vercellone et al. 



Giommi et al.ll2006l ). It also underwent major outbursts in 2007 (Ghisellini et al. 



2009a|) and 2008 (IVercellone et al. 



2010 



2006; 



Jorstad et al.ll2010l '). 



First observations of 3C 454.3 with the Fermi Large Area Telescope (LAT) began in July 
2008 during Fermi's commissioning p eriod, when the s ource was found at a high flux state with 
-^£;>ioo MeV — 10 X 10^^ ph cm^^ s^^ (JAbdo et al.ll2009l ). During this time, most observations were 



**School of Pure and Applied Natural Sciences, University of Kalmar, SE-391 82 Kalmar, Sweden 

* Corresponding author; cliarles.dermer@nrl.navy.inil 
Corresponding author: escande@cenbg.in2p3.fr 

* Corresponding author; lott@cenbg.in2p3.fr 
'Corresponding author; tanaka@astro.isas.jaxa.jp 



carried out in pointed mode with 3C 454.3 being close to the edge of the field of view (~ 55° off- 
axis), which did not allow for detailed spectral analysis during the brightest stage of the outburst. 
Observations in the decay stage, performed in survey mode, revealed a timescale less than 2 days 
for the flux to decline by a factor of 2. The spectrum showed a spectral break around 2 GeV 
with a spectral steepening from Fi = 2.3 to r2 = 3.5. Such a break has now been found to be a 



comm on feature in bright FSRQs and in some low-synchrotron peaked BL Lacs as well (lAbdo et al. 



2010dl ). Based on weekly light curves, a very moderate "harder when brighter" effect has also been 
observed, with the photon spectr al index obtained w ith a single power- law fit varying by less than 
0.3 for flux ratios varying by > 7 (JAbdo et~al]l2010dh . The source is hsted as IFGL J2253. 9+1608 



in the First LAT AGN catalog (jAbdo et al. 



2010 



3). 



The continuous monitoring by the Fermi LAT showed that the source activity faded contin- 
uously in early 2009 and then rose back up from June onwards. It underwent an exceptional 
outburst in November 2009 - January 2010 when it became the brightest 7-ray source in the sky 
for o ver a week, reaching a record daily flux level in the GeV band as seen by the LAT and AG- 



ILE (Striani et al. 



2009a 



b; Escande &: Tanaka 



show ed strong activity at optical frequencies ( Villata et al. 20091: Bonning et al. 



2009|), in the Swift /XRT and BAT bands (ISakamoto et al. 



200' 



i 



Striani et al. 



2OIOI'). At the same time it also 



200 



I 



2009; 



Krimm et al.l 



Sasada et al. 



20091 ^. and in the 



INTEGRAL/IBIS band JVercellone et al.ll2009bl '). iTavecchio et all (l20ld ') claimed 7-ray variability 
on timescales of a few hours from the LAT data and discussed its i mplic ations on the size and 
distance to the black hole of the emitting region, while iBonnoli et al.l (J2010l ) studied the single-day 
broad-band SEDs in the flare and in a lower state and found that a single-zone external Compton 
plus synchrotron self -Compton model adequately fltted the data. However, with a similar analysis, 
Pacciani et al.l (I2OIOI ) came to the conclusion that an additional component was required to flt the 
data. The source remained active afterwards with a slowly decaying flux around 2xl0~® ph cm~^ 
s~^ until early April 2010, when it brightened up again to a flux level of ~ 16x 10^^ ph cm^^ s^^, 
prompting the first Fermi-LAT target-of-opportunity (ToO) pointed observation beginning on 2010 
April 5 lasting for 200 ks. 

These two major events offer a unique opportunity to probe intraday variability and the as- 
sociated spectral changes in the 7-ray band, which is the focus of this paper. Particular attention 
is paid to the correlated spectral/fiux variations on different time scales, both in terms of spectral 
hardness and position of the spectral break, which have not been investigated in detail before. In 
Section 2, observations and analysis of data from 3C 454.3 from August 2009 through April 2010 is 
presented. Section 3 gives results of the analysis, and Section 4 provides interpretation. We sum- 
marize in Section 5. A flat ACDM cosmology with Hq = 71 km s~^ Mpc~^, $7^ = 0.27, fl\=0.73 
is used in this paper, implying a luminosity distance cIl = 1-69 x lO'^'^ cm to 3C 454.3. 



2. Observations and analysis 

The Fermi-hAT is a pair-conversion 7-ray telescope sensitive to photon energies greater than 
20 MeV. In its nomi nal scanning mode , it surveys the whole sky every 3 hours with a field of 



view of about 2.4 sr (jAtwood et al.ll2009l ). The data presented in this paper (restricted to the 100 
MeV-200 GeV range) were collected from 2009 August 27 to 2010 April 21 in survey mode, except 
for a 200 ks period starting on April 5 at 19:38 UT when the pointed mode was used, resulting 
in a gain of rate of accumulation of exposure by about a factor of 3.5 over survey mode. During 
pointed mode, the source direction was offset by 10° with respect to the LAT axis in order to limit 
adverse effects related to gaps in the detector that can affect on- axis photons. Over 3000 photons 
with E > 100 MeV were collected in pointed mode. To minimize systematics, only photons with 
energies greater than 100 MeV were considered in this analysis. In order to avoid contamination 
from Earth limb gamma-rays, a selection on events with zenith angle < 105° was applied. This 
analysis was performed with the standard analysis tool gtlike, part of the Fermi-LAT ScienceTools 
software package (version v9rl5p5). The P6_V3_DIFFUSE set of instrument response functions was 
used. This set includes a correction for the average reduction of the effective area due to pile-up 
effects as fewer photon events pass the rejection cuts. This correction is sufficient for integration 
times longer than a day. The residual energy and trigger-rate dependent acceptance variations 
present for shorter times as established with Vela and Galactic diffuse emission data (typically 
amounting to less than to 10% in the considered periods), have not been corrected for in this 
analysis. 

Photons were selected in a circular region of interest (ROI) 10° in radius, centered at the 
position of 3C 454.3. The isotropic background, including the sum of residual instrumental back- 
ground and extragalactic diffuse 7-ray background, was modeled by fitting this component at high 
galactic latitude (file provided with ScienceTools). The Galactic diffuse emission model version 
"glLiem_v02.fit," was used, with both flux and photon spectral index left free in the fit (the Galac- 
tic longitude and latitude of 3C 454.3 is 86.1° and —38.2°, respectively). All point sources lying 
within the ROI and a surrounding 5°-wide annulus with a flux greater than about 1% of the 
quiescent level of 3C 454.3 were modeled in the fit with single power-law distributions. 

Although the actual spectral shape is better reproduced by a broken power law, the spectral 
variations were investigated using the photon indices resulting from single power-law (PL) fits, 
which provided the best statistical accuracy. All light curves were produced using fluxes derived 
with power-law fits. 

Different analyses were performed by fitting the spectra with various models over the whole 
energy range covered by the LAT at E > 100 MeV: a broken power law N{E) = N(){E / Ehreak)~^' , 
with i = 1 if E < Ebreak and i = 2 \i E > Ehreak, a log parabola function 

iN{E) = Nq{E/Ep)-°'-^^°3'^^I^p'^ where Ep is fixed at 1 GeV) and a power law with exponential 
cutoff function {N{E) = Nq{E / Eq)~^ exp{— E / Ecutof f)) , and with a PL model over equally spaced 
logarithmic energy bins with the spectral index kept constant and equal to the value fitted over the 



whole range. 

In case of fits with broken power law (BPL) models, the break energy {Ebreak)^ which separates 
the photon energy ranges where different photon indices Fi (for E < Eij^eak) and r2 (for E > Ehreak) 
apply, could not be obtained directly from the fit because of convergence problems due to the non- 
smooth character of the BPL function at the break energy. It was instead computed from a 
log-likelihood profile fitting procedure, with a statistical uncertainty corresponding to a difference 
of — 2AL=1 in the log-likelihood function L with respect to its minimum. 

In order to minimize spurious correlations between flux and spectral index, the fluxes -Fe>Ei 
were also computed above the photon energy Ei where the correlation between integrated flux and 
index is minimal. This energy was derived from a power-law flt of the form 

N{E) = FE>Eo(r - 1) E~^/Ef+^ (with ^0=100 MeV). Since Fe>Ei = i^E>Eo(^l/^o)"^+^ min- 
imizing AFe>Ei /-^E>Ei with respect to Ei yields ln(£'i/£'o) = Cft/{F-e>'Eq C'rr), where Cpv and 
Crr are terms of the covariance matrix returned by the fit. In the 3C 454.3 case, Ei has been 
found to be 163 MeV over the ToO time range. The same value has been used for the other time 
periods as well. 

The estimated systematic uncertainty on the flux is 10% at 100 MeV, 5% at 500 MeV and 
20% at 10 GeV. The energy resolution is better than 10% over the range of measured Ebj-eafc- 



3. Results 

Figured] displays the daily light curves (red points) from MJD 55070 to 55307 (2009 August 27 
to 2010 April 21) for fluxes above 100 MeV. The periods showing the fastest flux variations during 
the December flare, with fluxes changing by more than a factor of 2 in amplitude, are enlarged in 
the insets, with the E > 100 MeV fluxes averaged over a daily, 6-hour, and 3-hour binning shown 
by the red, open blue and green data points, respectively. The error bars are statistical only. Three 
flares displaying a flux variation greater than a factor of 2 over less than a day (MJD 55167, 55170 
and 55195) have been studied extensively during this period. These flares exhibit v ery different 



proflle s and degrees of symmetry. The flrst two flares were fltted with the function (lAbdo et al. 



201031 ) 



F = 2Fo(e(*«-*)/^'- + e(*-*«)/^/)-^ + F,,,{t), (1) 

where Tr and Tf are the rising and falling times, Fq is the flare flux amplitude and Fiyg(i{t) is 
a (slowly varying) background flux, while for the third one a constant plateau flux was allowed 
between the rising and falling phases. The characteristic flare duration can be estimated as T^+Tf. 
We obtained Tr=0.37 d and T/=0.06 d for the MJD 55167 (Dec. 2) flare, and T^=0.07 d, T/-=0.26 
d for the MJD 55170 (Dec. 5) flare. The MJD 55195 flare, which is less prominent, shows a rise 
time of ~ 0.16 d. In the two early flares, the flux variation Fq was greater than 10^^ ph(> 100 
MeV) cm~^ s~^, and statistically significant factor-of-2 variations take place on timescales as short 
as 3 hours. 



In Figure [Tl the light curve of the July-August 2008 outburst (shifted by 511 days; see 



m 



Abdo et alj l2009l ) is shown for comparison. The resemblance of the two light curves is notable 



, although the estimated fluences are different by 40%. The brightest periods of the outbursts last 
for about 10 days, and are then followed by a long tail of fairly high activity upon which are super- 
imposed minor flares lasting for a few days. Although the outburst in April 2010 is longer, lasting 
for > 30 days, it exhibits similar patterns as the December 2009 outburst, where its maximum level 
is preceded by an intermediate elevated flux lasting for about 5 days. This feature could serve as 
an alert for an imminent surge in flux. 

The power density spectrum (PDS) of the 2 009 Nov. 5 - 2010 M arch 4 3h-light curve has been 



calculated in a similar way as that reported in lAbdo et al.l (|2010al ) for the 11-month light curve 



(3-day bin), normalized to fractional variance per frequency unit (/) and white- noise subtracted. 
This PDS indicates a power law power density, 1//° with index a = 1.50 it 0.06, i.e., intermediate 
between flickering (red noise) and shot noise (driven by Brownian processes). Such result confirms 
over a wider frequency range the value found with the first 11 months of data. In Figure [2] the 
first order structure function (SF), the PDS and the Morlet continuous wavelet transform (CWT) 
of the continuous, unprecedented-resolution, 7-ray light curve of 2009 Nov. - 2010 March are 
reported. A break around 6.5 days is suggested by the SF analysis, the power index slopes being 
a = 1.29 lb 0.10 between 3 hours and 6.5 days, and a = 1.64 it 0.10 between 6.5 and about 26 days, 
while at longer lags spurious drops due to the finite range become apparent. The PDS analysis 
confirms these values (a = 1.40±0.19 at high frequency and a = 1.56ib0.18 at low frequency). The 
temporal behavior of 3C 454.3 is therefore showing 1//^'^ universality from 3-hours to 11-months 
timescales. The 6.5-day break, consistently confirmed by the PSD and more clearly depicted by the 
SF, represents a steepening towards longer lags (flattening towards higher frequencies) and does not 
necessarily imply a local characteristic time scale, as it could simply be the point where two PDS 
components with different slopes are equally strong. Gamma-ray variability in 3C 454.3 can be seen 
as a short term realization of a stochastic mechanism, where structures that are resolved in shorter 
observations are simply averaged out in long observations, and where big outbursts correspond to 
statistical tails of the same process. 

The CWT in Figure [2] provides a local and detailed time series analysis, through the two- 
dimensional energy density function (modulus of the transform, filled color contour) computed 
using a Morlet waveform, providing the best tradeoff between localization and period/frequency 
resolution. For timescales below 1 day no local peaks are found, although some marginal features in 
this time range are found during the outburst state. The big outburst of December 2009 is, in fact, 
localized and decomposed in a chain of well-defined power CWT peaks. The 6.5 days timescale is 
confirmed by the major peak out of the cone of influence (localized at about MJD 55166, i.e, the 
onset of the outburst, MJD ~ 55166.2 - 55172). A second energetic peak in this period is found at 
about 2.5 days manifesting another dominating timescale during the outburst, while in the second 
period of the outburst two minor power peaks at about 19h and 1.3 days are also visible. Based 
on this CWT local analysis there is no evidence for structure on time scales shorter than about 12 



hours but shorter time scales close to the sampling scale cannot be ruled out. 

The 1 GeV daily light curve, shown in Fig. [3] top, closely resembles the 100 MeV light curve, 
hinting at little spectral variability. This behavior is confirmed by the very limited variation of the 
photon spectral index measured at E > Ei = 163 MeV, displayed in Figure [3] bottom by the daily 
average photon index (open blue symbols) as well as the weekly-averaged ones (solid black points). 
The near constancy of the spectrum is in accord with the re sults found from the July 2008 flare 



and the first 6-months of LAT data (lAbdo et al.ll2009l . l2010d|) 



The variation of the amplitude of the weekly photon indices is only Ar=0.35 (varying between 
2.35 and 2.7) during the period under consideration, but the variation is statistically significant. 
Comparing the weekly photon indices to their weighted average returns a reduced x^ =86.4/32, 
corresponding to a 10~ probability for a nonvariable signal. Particularly notable is a slight soft- 
ening of the spectrum during the periods of lowest fluxes, particularly during the periods of MJD 
55221-55235 and 55264-55278. There is also a suggestion that a progressive hardening over several 
weeks precedes a major outburst, but more such events will be needed to establish whether this 
behavior is typical. 

The correlation between daily spectral index and flux was further investigated by computing 
a correlation function for four segments of the November 2009 - April 2010 time range. In the first 
three segments (Dec 2009 flare, flare decay and inter flare region) we find evidence for a lag such 
that spectral index variations lead the flux by about 5 days. In contrast the correlation during the 
April 2010 flare is weak and if anything shows a lag in the opposite sense. Although the significance 
of the lag in these data is only about 2 cr we do see the same 5-day lag also in the data from the 
2008 flare. 

FigureHltop left (right) presents the weekly- (daily-) averaged power-law photon spectral index 
vs. flux above Ei. A weak "harder when brighter" effect shows up for weekly bins, but is almost 
washed out when considering daily bins. To make the trend clearer, photons were sorted in 3 
daily-flux bins {Fe^ <2.5, 2.5 < Fe-^ < 5, and 5 < Fe^ < 10, where Fe^ is the photon flux above Ei 
in units of 10^^ ph cm~^ s~^) and the analysis was repeated with the resulting photon files. The 
result, displayed as magenta points in the top right panel of Figure HI still exhibits a slight harder 
when brighter effect. The data point at Fe^ — 11 corresponds to the December 2"^^^ flare and is 
consistent with the trend observed at lower flux. The three other panels in Figure |4] show the photon 
spectral index vs. flux patterns for the three rapid flares, obtained either with a 6-hour binning 
(middle left and bottom left) or a 3- hour binning (middle right). The light curves at E > Ei (not 
above 100 MeV, as in Fig. 1) calculated at the time of the main flaring episodes are displayed in 
the insets of Figure [1] Despite the flux reaching the highest values for a non-GRB cosmic source, 
the statistical significance of these patterns is marginal, except for the MJD 55167 flare. Instead of 
a well-defined, universal pattern, a variety of patterns is found. The MJD 55167 fiare is associated 
with a clockwise pattern, with a fiux-doubling accompanied by an essentially constant (or weakly 
harder) photon spectral index. The reduced x^ for a constant fit of the photon spectral index 
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is 28.6/9. T he pattern is somewhat in dicative of a "hard-lag" effect possibly reflecting particle 



acceleration (jKirk Sz MastichiadisI Il999l ). The short MJD 55170 flare shows an indication for a 
counterclockwise loop {Xr — 6.0/7). The loop diagram for this flare on 3-hour timescales reveals 
the simplest ordered pattern, in this case a softening followed by a flux increase which then decays 
in flux and hardness. This could reflect the underlying timescale on which coherent variability 
takes place in 3C 454.3. Because this pattern does not recur in other flaring episodes, no strong 
conclusions can be made, however. The MJD 55195 flare shows some softening during the plateau 
(points 3-6) following the flux rise as expected from a cooling behavior, but the significance is low 
{Xr — 8.6/9). While it is difficult to draw any firm conclusion on acceleration and cooling from 
these patterns, the lack of strong spectral variability still provides clues to the underlying physical 
processes. 

Figure [5] shows the flux and photon spectral index as a function of time in the period around 
(blue symbols) and during (red symbols) the time of the target of opportunity when the Fermi 
LAT was in pointed mode (MJD 55291.82-55294.13). The binning is 6 hour and 3 hour for the 
survey and pointed modes respectively. As expected by the 3.5 fold increase in exposure per unit 
time during the ToO, the statistical accuracy in the measurement of both parameters improves 
significantly. Although in a high state, the source was unfortunately fairly steady during this 
period. No indication for variability more rapid than that observed during the giant outburst is 



found during the ToO period, as already noted bv lFoschini et al.l (J2010l ). The photon spectral index 



vs. flux above Ei measured in 3-hour time bins is plotted in Figure[6l The reduced x^ for a constant 
fit of the photon spectral index is 18.52/16, indicating that the data are consistent with a constant 
value. The correlation coefficient is 0.26. 

Figure [7] shows the integrated spectrum from Period 1, MJD 55121-55165, Period 2, MJD 
55166-55173 (week of the giant outburst). Period 3, MJD 55174-55262 and Period 4 MJD 55280- 
55300. The distributions have been fitted with a broken power law, a log parabola function and 
a power law with exponential cutoff function. The log parabola gives a worse fit than the broken 
power law and the power law with exponential cutoff functions, which are difficult to discriminate 
for these periods. The fitted parameters are given in Tabled] for the four periods. 

The variation of break energy (cutoff energy) with flux is displayed in the left (right) panel 
in Figure [8] at different observing periods. No strong evolution of either the break energy or the 
cutoff energy is found, but there is some indication of a slight hardening with flux. For a given flux, 
the position of the break energy is slightly different from that observed during the bright outburst 
in 2008, but all Eh rf.„k (and Ecutoff) are constant within a factor of ^ 2. For the same periods 



Striani et al.l (J2010l ) found fairly large spectral variability with the AGILE data, with a photon 
spectral index as low as 1.60±0.32 during big flares. This behavior is not conflrmed by the present 
analysis. 

The maximum photon energy found within the 95% containment radius from the location of 
3C 454.3 during the period MJD 55140-55261 was a photon with E = 20.7 GeV at MJD 55179.98, 
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when Fe>ioo MeV = 6 x 10^^ ph cm~^ s~^ and the variabihty time i^iaj. ^ 1 day. On December 
2, when the average of -Fe>ioo MeV over this day exceeded f« 20 x 10~^ ph cm~^ s~^, the highest 
energy was 8.5 GeV. As aheady noted, the flux variations during this flare were on timescales as 
short as a few hours or less. The energy of the most energetic photon during periods of rapid flux 
changes provides the strongest constraint on the 77 opacity, as discussed below. 



Discussion 



-2s-l 



from 3C 454.3 at z = 0.859 implies 



3.3 ± 0.15 X 10"^^ erg s'^. This is 



A photon flux of -F£;>ioo McV = 22 it 1 x 10^ ph cm 
an apparent isotropic 7-ray luminosity above 100 MeV of L^ 
~ 3 t imes larger than th e luminosity of the z = 1.839 blazar PKS 1502+106 during its August 2008 
flare (JAbdo et ahlboiOcI ). but stih lower than the luminosit y of PKS 1622-297 (^ 5.1 x 10"^^ erg s~i 
with the current cosmological model) during the 1995 flare (JMattox et al.lll997l ). The time-averaged 
7-ray luminosity (Tabled]) measured with the Fermi LAT is ~ 9 x 10^^ and ~ 1.6 x lO'^^ erg s~^ 
for periods 1 and 2, respectively. Estimates for t he black-hole mass in 3C 454.3 range from Mg ~ 1 
(JBonnoli et al.ll2010l ) to Mg fa 4 (JGu et al.ll200ll ). To be radiating below the Eddington luminosity 
of 1.26 xlO*^ (M9M0) erg s~^ for a lO^Mg Solar mass black hole implies that the high-energy 
radiation is beamed into a jet with opening angle 6j < 0.1 

The light curves from Fig. [1] show evidence for variability timescale of few hours and the 
model fitted to a very brief subflare at MJD 55170 implies a flux doubling timescale of (In 4) x T^. ~ 
2.3 h. However, if the bulk of the flare were characterized with such a short timescale, we would 
expect a very erratic light curve calculated with daily binning. Because the overall shape of the 
November/December 2009 flare is smooth, the dominating variability timescale is rather close 
to ivar ~ 1 d. The short subflares may reflect existence within or outside the main emitting 
zone of geometrical substructures as compact as R = (5ctvar/(l + -z) ~ 3.3 x 10^^525^var,2.3h cm ~ 

22(525 tvar,2.3h-Mgi?g. 

The minimum Doppler factor Smin < S = [T{1 — /3cos6)]~^, where 9 is the angle from the jet 
axis to the line of sight, is defined by the condition r^^(ei) = 1, and can be estimated to ~ 10% 
accuracy compared to results of more detailed numerical calculations through the expression 

1/6 

(2) 



aTdlil + zff.ei 



^^var'^eC- 



(Dondi &: Ghisellini 



1995 



Ackermann et al.ll2010l ). where /^ is the vF^ spectrum of 3C 454.3 mea- 
sured at frequency u = me(?t/h, and Ei = me(?ei is the maximum photon energy during the 
period in which /^ is measured. The vFy flux /^ in eq. ([2]) is evaluated at e = e = 25^/(1 + z)'^ei 
from the pair-production threshold condition. Writing /^ = 10~^''/_io 6rg cm~^ s~^ in eq. ^ gives 
5min ~ 13 [/_io^i(10 GeV)/tvar,d]^^'^ and E{ keV) ^ 3.4((5min/15)V^i(10 GeV). Contemporaneous 
Swift XRT observa tions give the i/F,y fl ux between 0.2 and 10 keV, with i^Fj,(4 keV) ?s 6 x 10^^^ 



erg cm ^ s ^ (e.g.. iBonnoli et al.ll2010l ). so that 5mm ~ 15. 
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The value of the Doppler factor can also be constrained from the argument that the energy 
flux from the Synchrotron Self-Compton (SSC) component should be below that observed at soft 
X-rays, in terms of isotropic luminosities Lssc ^ -^x- SSC luminosity is related to the synchrotron 
(SYN) luminosity via 1/ssc/-^SYN ~ '^syn/''^b ^^ ^^^ Thomson limit, where u'^yn ~ -^syn/(^'^-^^'^) 
is the energy density of synchrotron radiation and Ug = S'^/Stt is the energy density of magnetic 
field. Noting that Lorentz transformation of bolometric luminosity is Lsyn 



(5^LgYN) we obtain 



8'> 



2{i + zy^YN 



(3) 



This can be combined with a constraint from the external-Compton (EC) and synchrotron lumi- 
nosities ratio, 



L 



SYN 



s-2 ItEXT 



"B 



(4) 



(JDermeiill995l ). where uext is the energy density of external radiation. Eliminating from these 
equations the magnetic field value, we obtain 

-1/8 

(5) 



-1/8 



-1/8 



^ ll.O X i^gYN,48 ^ -^^£0,49 ^ -'^X,47 ^ ''var.d ^ 



UEXT 



0.015 erg cm 



On the other hand, we can eliminate the Doppler factor and obtain a lower limit on the magnetic 
field 



-5/8 



-3/8 



1/4 



B > 3.4 G X I'gYN,48 ^ ^ECA9 ^ -^X,47 ^ *var,d 



^EXT 



3/8 



_.. ; • (6) 

0.015 erg cm •^y 

Order-of- magnitude r eference values of iso tropic luminosities -^syn, Lx and Leg have been de- 
duced from Fig. 4 in iBonnoli et al.l (J2010l ). If the source region is located within the broad- line 
region (BLR), the energy density of external radiation is uext ~ -^BLR/(47rrgLRc), where Lblr is 
the luminosity of the broad emission lines and tblr is the characteristic radius of the BLR. We 
approximate -Lblr with GALEX measurement of Lya line, Li^y a ~ 2 x 10^ ^ erg s~^ and adopt 



^RT.F! ^ 6 X 10^ cm based on reverberation scaling relations for CIV line ( Bonnoli et al. 



Kaspi et al.ll2007l ^. 



201C 



The analysis of radio observations made at a different epoch gives 6 = 24.6 it 4.5, bulk Lorentz 
fa ctor r = 15.6 ifc 2. 2, and observing angle = 1.3° it 1.2° obtained from super luminal observations 



(jjorstad et al.ll2005l ) . Both our constraints on the Doppler factor are consistent with this result and 
would become even tighter, if we adopt a shorter variability timescale. Hence, we take 625 ~ r2o ~ 1 
(r = 20r2o)- The constraint on the value of the magnetic field B' is a few times the equipartition 
magnetic field B',^{G) « 2 (/.loAs)'/^ (4'^''^iT*w,d)> ^here /sYN,peak = 10-i%io erg cm-^ 
s~^ is the observed synchrotron peak flux, z^sYN.peak = 10^'^ 1^13 Hz is the synchrotron peak frequency 
and A is the prod uct of the bolomet ric correction factor and ratio of energy contents in protons 
and electrons (e.g.. iFinke et al.ll2008l ). 



For a conical geometry of opening angle 9j > R/r ~ l/F, the location of the emitting region 
for the December flare is constrained to be at distance r < 2cr^ivar/(l + z) !^ 0.2rf5ivar,d pc; i.e.. 



13 



towards the outer parts of the BLR. A narrow jet with opening angle < 0.2/r , places the emitting 
regio n further out, and could be consistent with a fragmentation of outflow (JMarscher &: Jorstad 
20ld ). The most rapidly varying flaring episodes suggest locatio ns at the sub-pc sc ale, contrary to 
conc lusions from coherent optical polarization changes in 3C 279 (jAbdo et al.ll2010d ) and PKS 1510- 
089 ([Marscher et al.ll20ld ) over timescales of tens of days, which suggest that the emitting regions 
are several pc from the black hole. Larger distances would be inferred if L ^ 15 or 6j <^ l/L. Higher 
r factors would conflict with the expected source density of th e parent population of FSRQs, as we 



discu ss later. Narrow jets are suggested by radio observations (jPushkarev et al 



200S 



Jorstad et al, 



20051) and could result frorn recollimation shocks at the pc scale (e.g., iNalewajko &: Sikorall2009l : 
Bromberg &: LevinsonI l2009l ) that reduce the characteristic size of the emission region. A multi- 
wavelength radio, mm, optical, and X-ray campaign covering observations of 3C 454.3 between 
2005 and 2008 shows coherent optical and mm flux and polarization episodes consistent with the 
bright est events taking pl ace at the mm core located ~ pc beyond the acceleration and collimation 
zone (| Jorstad et al.ll2010l ). Without strong collimation or recollim ation, the short obser ved vari- 



ability timescale implies a very small value 0,- ~ 2 x 10 ^tvar.d/'^pc ( Tavecchio et al.ll20ld ) with an 



associated low radiative efficiency. Alternately, flaring episodes with short variability times might 
take place within the BLR, whereas the more slowly varying emissions could be radiated by jet 
plasma at larger distances. Spectral variations due to the different target photon sources might be 
concealed by mixing from components made at small and large radii. "Residual" collisions at large 
radii might also m ake a slowly va.rying underlying emission components, as similarly inferred for 
GRB observations (JLi fc WaxmaiJl2008l ). 



The large changes in flux exhibited by the 7-ray light curves in Fig. [T] could be due to several 
different factors, including a changing mean magnetic field, electron number, Doppler factor, target 
photon density and spectrum as the jet moves outward, or some combination of these factors. Any 
such explanation must also account for the moderate spectral changes, including the near constancy 
of -Ebreak (Fig. [8|), while the flux changes by over an order of magnitude. 

One possibility is related to scattering of a target photon field in the Klein-Nishina regime. 
Compton scattering takes place in the Thomson limit when the energy of the photon to be scat- 
tered is (in mgC^ units) e" < 1/4 in the electron rest-frame. Since the radiation from ultrarela- 
tivistic electrons is seen only when they are moving very close to the line of sight, that energy 
for 7 ^ r is e" ~ 'jSe^,, where e* is the average energy of the target photons in the lab frame. 
Hence, the scattering is in the Thomson regime provided the upscattered photon energies are 
ec ^ {i/3)6^'y^ej{l + z ) < l/\12ejl + z)], i.e. for £;c(GeV) < U/E^eV), independent of the 



Doppler factor (compare iGeorganopoulos et al.ll2001 



If the break energy observed in 3C 454.3 at ~ 2 GeV is due to the transition to scattering 
in the KN regime, then the underlying target photon energy i?^ ~ 6 eV is close to the energy of 
the Ly a photon at 10.2 eV. We have tested this possibility by comparing the Compton-scattered 
spectrum from a power-law electron distribution with a monochromatic Ly a photon source with 
the Fermi LAT data. The spectrum is independent of 6, though can depend on maximum and 
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niinimum electron Lorentz factors. As shown in Fig. [9l the spectrum is too hard to fit the data, 



and t reatment of KN effects on coohng (JDernier &: Atovanll2002l : lModerski et al.ll2005l : ISikora et al 



20091 ) or the addition of other soft photon sources will reduce the sharpness of the break. The 
difficulty to fit the sh arp spectral break with a single power-law electron distribution is in accord 
with the conclusion of lAbdo et al.l (|2009l ) that this break refl ects a complex electron spectrum. To 
obtain a good spectral fit to the 7-ray spectrum of 3C 454.3, iFinke &: Dermed (J2010l ) use a broken 
power-law electron distribution and a multi-component Compton-scattering model. The robustness 
of Efyf^ak in this model is due to similar, ~ r~'^ dependence of accretion disk and external radiation 
energy density within the BLR. 



Although lAbdo et al.l (J2009) excluded the possibility of explaining the spectral break with 



photon-photon pair absorption, iPoutanen fc SternI (I2OIOI ) proposed that significant intrinsic ab- 
sorption can be provided by He II Lyman-a line at 40.8 eV and continuum at 54.4 eV. However 
their model implies the location of the gamma-ray production is in the inner edge of the BLR 
and inefficient scattering in the KN regime is predicted to lead to significant hardening of syn- 
chrotron spectra in IR/optical bands. Further studies are required to verify whether such effect is 
in contradiction with multiwavelength spectra created within a one-zone model. 



15 and the jet opening angle 6j ~ l/P, then the beaming factor for a two-sided top-hat 
: l/2r^ = l/(450rfc;) . Within the framework of the unification hypothesis for radio 



If r 

jet is fb 

galaxies (lUrry fc Padovanilll995l ). FRII radio galaxies are the parent popula tion of FSRQs. Th e 

space density of FRII radio galaxies at 0.8 < z <1.5 is p = 3 x 10^'' Mpc~^ (JGendre et al.ll20ld ). 



At redshift unity, the volume of the universe 



A7r/3R% ^ 



3 X 10^^ Mpc^. Therefore, if ah FRII 



radio galaxies were misaligned FSRQs similar to 3C 454.3, then the implied number A^ of FSRQs 
is N < fhX3 X lO"'^ Mpc"3 X 3 X 10^^ Mpc^ ^ 200/r?^,. T his number is similar to the number, 
~ 100, of 7-ray blazars at redshift unity (jAbdo et al.ll2010bl ). so is in accord with the unification 



scenario. This simple argument suggests that unless 9j ^ l/F, in which case the total jet energy 
of an individual FSRQ blazar would be correspondingly larger, the typical bulk Lorentz factor F 
for FSRQs would not be > 15. 



5. Summary 

The correlated spectral and temporal properties of 3C 454.3 during two very strong fiaring 
episodes, during which the source was the brightest object in the 7-ray sky, have been studied. The 
7-ray spectrum of 3C 454.3 shows a significant spectral break between ~ 2 - 3 GeV that is very 
weakly dependent on the flux state, even when the flux changes by an order of magnitude. The 
light curves during periods of enhanced activity in July - August 2008 and December 2010 show 
strong resemblance, with a flux plateau of a few days preceding the major flare. The spectral index 
measured on a daily basis shows a very moderate "harder when brighter" effect for fluxes measured 
at E > El, where Ei = 163 MeV is chosen to minimize spurious correlations resulting from the 
choice of the low-energy bound on the range. However, the weekly spectral index displays a more 
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significant variation, but not exceeding Ar=0.35 when the flux varies by more than a factor of 
10. Indication for a gradual hardening preceding a major flare and extending for several weeks is 
found. No recurring pattern in the photon spectral index/flux plane measured at E > Ei during 
the course of flaring episodes has been identified. Flux variations of a factor of 2 have been observed 
over time scales as short as a few hours, though only weak variability was observed during the time 
of the Target of Opportunity pointing of the Fermi Telescope towards 3C 454.3. 

From 77 opacity constraints, we derive a minimum Doppler factor 6min — 15 from the fiux, 
variability time, and highest energy photon measurements. T his value is in accord with independent 



measurements of 5 from superluminal motion observations ([Jorstad et al.ll2005l ). The behavior of 
the break energy has also been investigated. Spectral softening due to the onset of Klein-Nishina 
effects in Compton scattering, which is independent of Doppler factor for a power-law electron 
distribution, was considered as the origin of the 7-ray spectrum. The magnitude of the softening 
for such a model was found to be inadequate to fit the spectrum without introducing intrinsic 
breaks in the underlying electron spectrum. An independent estimate of 5 and B also gives values 
of 5 > 15; also B > several G. The short fiaring times suggest an origin of the emission region 
within the BLR, but uncertainties in parameters are compatible with the production of emission 
from 3C 454.3 on the pc scale. 
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Fig. 1.— The light curve of the flux of 3C 454.3 in the 100 MeV ~ 200 GeV band (red). The hght 
curve of the July - August 2008 flare, shifted by 511 days, is shown for comparison (magenta). The 
insets show blow-ups of the two periods when the largest relative flux increases took place. The 
red, blue, and green data points in the insets correspond to daily, 6 hr, and 3 hr averaged fluxes, 
respectively. The fit results discussed in the text are displayed as solid curves. 
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Fig. 2. — Top panel: structure function of the 3h-bin flux light curve for the period 2009 Nov. 5 - 
2010 Mar. 4 and corresponding power density spectrum (inset). Bottom panel: plane contour plot 
of the continuous Morlet wavelet transform power density for the same light curve. Thick black 
contours are the 90% confidence levels of true signal features against white/red noise background, 
and cross-hatched regions represent the "cone of influence", where spurious edge effects become 
important. 
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Fig. 3. — Top: Light curve of the flux in the 1 GeV - 200 GeV band. Bottom: Variation of the 
daily (blue points) and weekly (black points) photon spectral index derived from a power-law fit. 
The black line depicts the mean weekly spectral index. 
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Fig. 4. — Photon spectral index vs. flux measured on a weekly (upper left) and daily (upper 
right) basis are shown by blue points. The averages over periods with daily fluxes Fei < 2.5, 
2.5 < Fei < 5, 5 < Fe^ < 10, where Fe^ is the photon flux above Ei in units of 10~^ ph cm^^ 
s^^, are shown by the magenta data points in the upper right panel. Middle left, middle right, and 
bottom panels show photon index vs. flux for the MJD 55167, MJD 55170 and MJD 55195 flares, 
respectively. The labels refer to the different times of the light curve as shown in the insets. 
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Fig. 5. — Flux (filled data points; left-hand axis) and photon index (open data points; right-hand 
axis) as a function of time in the period surrounding the ToO pointing. Data collected in survey 
mode (6-hour binning) are in blue, those collected in pointed mode are in red (3-hour binning) 
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Fig. 6. — Photon spectral index vs flux on 3-hour timescales during the time of the ToO pointing. 



23 




Energy (MeV) 




Energy (MeV) 



Energy (MeV) 



Fig. 7. — uFy distributions of the Fermi LAT data for different time periods, along with the 
fitted broken power law (solid), log parabola (dashed) and power law+exponential cutoff (dotted) 
functions. 
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Fig. 8. — Left: E'broak vs flux for different time periods (given in MJD). Right: E'cutoff vs flux for 
the same time periods. 
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Fig. 9. — Model for the 7-ray spectrum of 3C 454.3 when a power-law electron distribution 
Compton scatters Ly a photons. Best fit electron number index p = 3.17. Model is insensitive to 
values of lower and upper comoving electron Lorentz factors 7mm and "fmax provided 'fmin ^ 10^ 
and ^max ^ 10^. The Klein-Nishina softening from a power-law electron distribution gives a poor 
fit to the data. 
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period 


Ti 


r2 


'^break 

(MeV) 


a 


P 


r 


^cutoff 

(MeV) 


1 


2.31± 0.02 


3.19±0.12 


1810ti° 


2.52±0.03 


0.11±0.02 


2.29 ±0.03 


5200±900 


2 


2.33± 0.03 


3.29±0.21 


2750+t- 


2.48±0.03 


0.09±0.02 


2.25±0.04 


7000±1700 


3 


2.32± 0.02 


3.12±0.08 


1630t\ll 


2.54±0.02 


0.12±0.01 


2.24±0.03 


5100±620 


4 


2.38± 0.02 


3.23±0.14 


238011° 


2.55±0.03 


0.10±0.02 


2.30±0.03 


7000±1300 



Table 1: Parameters of the functions fitted to the spectra for the different periods considered in 
Figure El 



